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AN ATLAS OF SYNTHETIC LINE PROFILES OF PLANETARY NEBULAE 
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RESUMEN 
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We have constructed a grid of photoionization models of spherical, elliptical and 
bipolar planetary nebulae. Assuming different velocity fields, we have computed 
line profiles corresponding to different orientations, slit sizes and positions. The 
atlas is meant both for didactic purposes and for the interpretation of data on real 
nebulae. As an application, we have shown that line profiles are often degenerate, 
and that recovering the geometry and velocity field from observations requires lines 
from ions with different masses and different ionization potentials. We have also 
shown that the empirical way to measure mass-weighted expansion velocities from 
observed line widths is reasonably accurate if considering the HWHM. For distant 
nebulae, entirely covered by the slit, the unknown geometry and orientation do 
not alter the measured velocities statistically. The atlas is freely accessible from 
internet. The Cloudy _3D suite and the associated VISNEB tool are available on 
request. 

ABSTRACT 

We have constructed a grid of photoionization models of spherical, elliptical 
and bipolar planetary nebulae. Assuming different velocity fields, we have computed 
line profiles corresponding to different orientations, slit sizes and positions. The 
atlas is meant both for didactic purposes and for the interpretation of data on real 
nebulae. As an application, we have shown that line profiles are often degenerate, 
and that recovering the geometry and velocity field from observations requires lines 
from ions with different masses and different ionization potentials. We have also 
shown that the empirical way to measure mass-weighted expansion velocities from 
observed line widths is reasonably accurate if considering the HWHM. For distant 
nebulae, entirely covered by the slit, the unknown geometry and orientation do 
not alter the measured velocities statistically. The atlas is freely accessible from 
internet. The Cloudy _3D suite and the associated VISNEB tool are available on 
request. 
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1. INTRODUCTION 

As is known, nebular line profiles bear informa- 
tion on the motions of the emitting gas. The study 
of line profiles thus allows one to measure the expan- 
sion velocity of nebulae and even to determine their 
internal velocity field. The velocity field is an im- 
portant key for understanding the dynamics and the 
genesis of planetary nebulae. The expansion velocity 
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allows one to derive the expansion age, and from this, 
such parameters a s the masses of the central stars 
(iGornv et al.l 119971 ). However, the determination of 
these characteristics is a very difficult task in gen- 
eral, since only the velocity field perpendicular to the 
plane of the sky can be observed, through the study 
of line profiles. Even if one could, for a series of ions 
with different ionization potentials, obtain line pro- 
files in every pixel of a high-resolution image of a PN, 
one would still need to make assumptions about the 
projection of the velocity vectors in the plane of the 
sky, in order to reach a full description of the veloc- 
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ity field. Fortunately, most PNe exhibit some kind 
of overall regularity in their morphology, and some 
degree of symmetry (apparently circular, ellipsoidal, 
bipolar...). So, if an image of the nebula is available, 
it is possible, with the help of a few assumptions, to 
determine the entire nebular geometry and velocity 
field. Most of the time, however, only limited data 
are available: line profiles obtained through one slit, 
or even, if the object has very small angular dimen- 
sions, line profiles corresponding to the emission of 
the entire object. If the object is faint, line profiles 
can be obtained only for the brightest lines, the ex- 
treme case being when they are available only for 
[O III] 5007A. Finally, morphological information on 
the object may be completely lacking, as is the case 
for PNe in distant galaxies. In such cases, what is 
the value of the information that one can derive from 
line profiles? 

In practice, the analysis of line profiles of plan- 
etary nebulae has followed several paths: i) A 
simple determination of the expansion velocities 
from observed line widths or from line sp litting 
(most of the references quoted in IWeinberger|[l98 9) : 
i i) Fitting of line profiles with spherical models 
([Gesicki et al.ll2QQ3f ): iii) Modeling e mission line pro- 
files and Positio n- Velocity diagrams ("Sabbadin^igS^; 
IWeedmanl 1 1968) and more recently Meaburn et al. 



(2005) and ISteffen fc Loped (|2006D : ISteffen et all 
J2OO7) which their powerful tool SHAPE; iv) Fitting 



of line profiles and other nebul ar characteristics with 
a 3D photoionization model (iMorisset et al. I I2OOOI : 
iMonteiro etlo] 1200 0": 'SabbadinlllQSj; v) Tomogra- 
phy ([Sabbadin et al. 2006, and references therein). 

The question of the robustness of the results is 
difficult and is almost never addressed in such stud- 
ies. The atlas of line profiles that we present in this 
paper is meant as a tool to visualize the variety of 
possible line profiles, depending on the morphology 
of the nebula, the velocity field and the observing 
conditions. The atlas is to be installed in a virtual 
observatory environment, that will collect in a large 
data base the models and line profiles com puted by 
any u ser with the tools we have developed (|Morissetl 
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The description of the atlas is given in Sec. [2l In 
the second part of this paper, we use the atlas, in 
its present stat, for two applications: a search for 
models reproducing a given profile (Sec. [3|) and an 
estimation of the robustness of the determination of 
expansion velocities based on profiles (Sec. [4]). 

2. DESCRIPTION OF THE A TLAS 

20061 



We use Cloudy_3D (jMorisset 



here- 
after C3D), which is a pseudo-3D photoionization 
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Fig. 1. Figure defining the geometry and the parameters 
a and h. 

code, similar t o the one presented in detail in 
Morisset et al. ( 2005L but based on the ID pho- 
toionization code Cloudy (Ferland eLalJ [l998) in- 
stead of NEBU. The computations for this paper 
were based on version 07.02.01 of Cloudy. The size 
of the cube in which the interpolation of the Cloudy 
ID models is done is set to 150^. We run ten Cloudy 
ID models in different directions. Once a 3D pho- 
toionization model has been obtained, a velocity law 
is applied to compute line profiles using the section 
VIS NEB of the C3D pa ckage, in the sa me way as de- 
scribed in lMorisset fc Stasihska ( 20061 ). The profiles 
are computed on a 55 wavelengths grid covering the 
range [-60, +60] km.s~^. 

The fine structures of the HI, Hel and Hell lines 
are taken into account by splitting the lines into two 
components. The relative intensities, velocity shifts 
and FWHM corrections for the two components are 
described in Tab . [H The values for Hell 4686A are 
adapted from'Clegg et al.l (|l999l ) and P. Storey (pri- 
vate communication). 

2.1. The parameter space and the nomenclature of 
the models 

All the models in the present version of the at- 
las have been run by changing the parameters which 
have the largest impact on line profiles: i.e. morphol- 
ogy and velocity law. The models are viewed with 
different orientations and different slit sizes and loca- 
tions, all of which strongly affecting the observed line 
profiles. On the other hand, parameters that play 
only a secondary or even negligible role have been 
kept unchanged. This concerns the effective temper- 
ature of the central star (chosen to be T*=10^K, so 
as to produce low and high ionization ions), the stel- 
lar luminosity (set to 10^ L0), the chemical composi- 
tion (chosen to be the Cloudy default), the absolute 
value of the density. In this version of the atlas all 
the computed models are ionization-bounded. 

At each point of the nebula with r > r^n, the 
hydrogen density is determined by the following law: 
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TABLE 1 

PARAMETERS FOR THE RECOMBINATION FINE STRUCTURE LINES. THE GLOBAL SHAPE IS A 
SUM OF TWO GAUSSIAN FUNCTIONS. II AND 12 ARE THE RELATIVE INTENSITIES OF THE 

TWO COMPONENTS, SHIFTl AND SHIFT2 THE POSITIONS OF THE CENTER OF EACH 
GAUSSIAN RELATIVE TO THE CENTROID, SIGMAl AND SIGMA2 THE FACTOR APPLIED TO 
THE NATURAL+THERMAL BROADENING TO TAKE INTO ACCOUNT THE VARIOUS 
SUBCOMPONENTS IN EACH GAUSSIAN (ONLY FOR HeU 4686 A). 
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Fig. 2. Figure defining the x, y, and z axes and showing 
which is the symmetry axis. 



where Ar(i9) = 1.6 x Rstr{0) with Rstr{0) the 
Stromgren radius in the direction defined by the po- 
lar angle ^ = in the equatorial direction). The 
first term of Eq. ([T|) is the inner equatorial density 
fixed to 2000 cm~^, the second term refers to 
the angular dependency of the inner density using 
a parameter p, and the third term is the gaussian 
dependency of the radial density. The inner shape 
of the nebula defined by r^n (0) is determined by two 
parameters a and h shown in Fig. [TJ where a is the 
eccentricity of the ellipsoids and h is the distance 
between the center of the nebula and the centers of 
the two lobes (a value of 0.0 for h means a purely 
ellipsoidal or spherical nebula). 

The adopted nomenclature of the models is the 
following: 

Irst three digits - generic name of the atlas, in 
the present version : PN3. 

2nd two digits - overall geometry of the inner 
cavity (defining ViniO)): SP: sphere (a=l, /i=0); 
EL: ellipse (a=1.5, h={))] BS: bipolar from two 
spheres (a=l, /i=0.5); BE: bipolar from two ellipses 
(a=1.5, /i=0.5); BL: bipolar from two elongated el- 
lipses (a=2.2, /i=0.5). 

3rd digit - parameter p in Eq. ([1]). 

4th digit - size of the central cavity. S: small (i.e. 
corresponding to f=20% of the Stromgren radius in 
each direction); L: large (i.e. f=80%). 

5th digit - velocity law. B: ballistic flow (i.e. v ex 
r); C: constant velocity. 

6h digit - turbulence. 1: turbulence of 10 km.s~^; 
0: no turbulence. 

7th digit - orientation. F: face on; P: pole on; T: 
tilted by 45 degrees (See Fig. [2]). 

For example, the model PN3_BE_1_L_B_0_T cor- 
responds to a bipolar nebula obtained from 2 ellipses, 
with a decrease of the density following the inverse 
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of the inner radius of the nebula, with a large cen- 
tral cavity, expanding with a ballistic velocity law 
without turbulence, seen tilted by 45 degrees. 

In its present version, the atlas is thus composed 
of 26 photoionization models corresponding to 5 ge- 
ometries, 3 angular density laws and 2 cavity sizes 
(the angular density law does not change the spher- 
ical models). For each model, four velocity fields are 
considered. So, in total, there are 104 virtual nebu- 
lae, each of which can be observed from 3 different 
directions (which are identical in the case of spheres, 
leading to 308 pages for the atlas). 

Figure [3] shows the images of all the models in the 
atlas, with their various orientations, co-adding the 
intensities of the H/3, [Nil] 6584A and [Olll] 5007A 
lines. This "model chart" allows a quick orientation 
in the atlas. 

The virtual observations of these models are 
made through 4 different kinds of apertures: a small 
square aperture, a long vertical slit, a long horizontal 
slit, and an aperture covering the entire object. We 
also compute the case where the aperture is not cen- 
tered on the object, leading to a total of 7 apertures 
(the big aperture is not affected by the off-center 
shift). In all cases, the line profiles were computed 
with a seeing of 1.5''. The seeing needs to be consid- 
ered during the computation of the intensities pass- 
ing through the selected apertures. On the other 
hand, we chose to not apply any instrumental pro- 
file, as this effect can be taken into account by any 
user downloading the data from the server. 

All the nebulae are set to be at 3 kpc from the 
observer. On the other hand, the geometrical size 
of the cube in which the nebula is constructed by 
C3D is given by the maximum extension of the neb- 
ula, which is controlled by the outer radius in the 
polar direction. This value is changing from one 
model to the other. In the case of elongated neb- 
ulae viewed pole-on, we basically see the equatorial 
torus, which appears to be smaller than in the case 
of spherical nebulae (see Fig. [3]). Taking the same 
distance for all the nebulae leads to changes in the ra- 
tio pixel/angular size from one model to the other in 
such a way that the relative size of the slits compared 
to the size of the equatorial torus is kept constant. 
Each virtual observation delivers profiles in 5 lines 
(namely: R(3, He I 5876A, Hell 4686A, [Nil] 6584A 
and [O III] 5007A). A total of 10780 different profiles 
are thus available (resulting from the computation 
of about two billions individual spectra!). 

2.2. Description of one page of the atlas 
Each page of the atlas (an example is given in 
Fig. [4]) corresponds to a different geometry, veloc- 
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Fig. 3. The model chart: images of all the models in a 
the sum of the light in the H/3, [N II] 6584A and [O III] 
5007A lines. 
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ity law and orientation. The first row describes the 
model parameters. Panel 1 shows the distribution 
of riH as a function of r/vmax {'^max being the half 
size of the cube, corresponding to the polar exten- 
sion) along the small axis (continuous line) and along 
the large axis (dashed line). Panel 2 shows, with 
the same conventions, the distribution of the elec- 
tron temperature Tg as computed by C3D. Panel 3 
illustrates the assumed velocity field, again with the 
same conventions. Panels 4 and 5 give additional 
information on the model definition. 

The panels in the second row represent 
monochromatic images (for the orientation defined 
in Panel 4) in the lines He I 5876A, He II 4686A, 
[N II] 6584A and [O III] 5007A. The positions of the 
slits are also indicated, in solid lines for centrally 
positioned slits, in dashed lines for non central slits. 

In all the following rows, the emission lines con- 
sidered are the same as in the second row, and they 
are ordered in the same way, from left to right. 

The panels in the third row show the surface 
brightness distribution in those lines. The contin- 
uous curve corresponds to the small axis, the dotted 
curve to the large axis. 

The next four rows show the line profiles, with 
the bold curve corresponding to a centrally posi- 
tioned aperture, the thin one corresponding to a 
non central aperture. The row number 4 (panels 16 
through 20) corresponds to a small slit of 3''x3''. 
The row number 5 (panels 21 through 25) corre- 
sponds to a horizontal slit of 3^' width, while the 
row number 6 (panels 26 through 30) corresponds 
to a vertical slit of same width. Finally, the row 
number 7 (panels 31 through 35) corresponds to an 
aperture covering the entire nebula. As expected, 
the line profiles show more complex structures if the 
slit is small. 

Rows 8 and 9 show position velocity diagrams 
(or echellograms) obtained through a vertical slit 
centered on the nebula (panels 36 through 40) and 
through an horizontal slit centered on the nebula 
(panels 41 through 45). 

Finally, row 10 (panels 46 through 50) shows the 
channel maps of the gas with null radial velocity. 

In rows 4 through 7, the horizontal segments indi- 
cate the values of the expansion velocities measured 
by classical methods (for centered slits only). Red: 
half peak-to-peak velocity. Blue: expansion veloc- 
ity measured from the half width at half maximum 
(HWHM). Green: expansion velocity from the half 
wi dth at one te n th ma ximum (HWIOM), as given 
by iDopita et al.l ( 1985 ). The black segment gives 
the mass-weighted expansion velocity perpendicular 



to the plane of the sky (actually the quadratic sum 
of the expansion velocity, the thermal velocity and 
the turbulent velocity, weighted by 11^+). For each 
segment the corresponding velocity is also indicated 
on the right side of the plot, with the mass- weighted 
velocity using bigger character sizes. 

2.3. Browsing the atlas 

The entire atlas (images, figures sim- 
ilar to Figs, m and [5] and emission line 
profiles in ascii format) is available at 
http : //132 . 248 . 1 . 102/Atlas_prof il es/, 

The atlas exhibits a large variety of line pro- 
files. In general, the profiles showing the most com- 
plex structures are those obtained through long slits. 
This is because long slits probe kinetically differ- 
ent regions, while the smoothing effect is minimized. 
Non central slits always produce asymmetric line 
profiles. In the case of the He I 5876A and Hell 
4686A lines, the presence of multiple fine structure 
components also induces asymmetric profiles, even 
for symmetric nebula and centered slits. The most 
complex line profiles generally correspond to pole- 
on nebulae. Some of those profiles are re miniscent 
of those observed bv iGesicki et al.l ( 20031 ). As ex- 
pected, double peak profiles principally occur if ob- 
servations are made with small slits while large slits 
give single peak profiles. In general, ballistic expan- 
sion tends to produce triangular profiles (except for 
[N II] 6584A which is emitted in the outskirts of the 
nebula and may produce a square profile) while uni- 
form expansion associated with turbulence tends to 
produce bell-shaped profiles. 

However, similar profiles may correspond to very 
different models, as will be better shown in the next 
section. 

3. WHICH MODELS CAN ACCOUNT FOR A 
GIVEN LINE PROFILE? 

So far, in the literature, most of the analy- 
sis of line profiles has be en either purely empirical 
(e.g. iMedina et al.l 2006) or by fitt ing to a simple, 
spherical model (^Gesicki et al. I I2OO3I . and references 
therein). However, a large fraction of PNe for which 
images are available are not spherical. 

One can then wonder how well the velocity field 
inferred from a model fitting with a spherical neb- 
ula represents the true velocity field and whether 
other solutions, implying non spherical nebulae, 
may fit the o bserv ed data as well. For example, 
Gesicki et al.l (|2003l ) have invoked the presence of a 
turbulent velocity field in order to fit their spherical 
models to the data. 
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Velocity field 
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Model: PN3_BE_1 _L_B_0_T 
Fri Jan 1 1 20:04:12 CST 2008 
T* = 10^ K, L* = 10^ Lsun 
R max = 1.20e+18 cm 
Solar Abundances 
Orientation: 45. 0. 0. 
Projection on z 



Geometry: 

f = 0.80 h = 0.50 a = 1.50 p=1.00 
S=1.60 Nh°=2000 Log Qo=47.86 
Velocity field: 

V=0.0 + 40.0 X R (Vt,,b=0.5) 
Seeing = 1.5", I.R. =0.0km/s 
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Fig. 4. One page of the atlas, see text for detailed description. 
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Our atlas allows us to do the following experi- 
ment. We chose one line profile from the atlas, for 
example obtained through a slit covering the entire 
object and we search which virtual observations in 
the atlas give similar line profiles. For this, we define 
the difference between two profiles Pi (A) and P2 (A) 
as Di,2 = EW{X) X (Pi(A) -P2(A))VE^(A), for 
example giving the same weight W(A) = 1 to all the 
wavelengths (as apparently done by Gesicki). 

The result of such an experiment is shown in 
Fig. m Each row corresponds to a different emis- 
sion line (see caption). The right panels show the 
profiles (in the same line) that differ by Di^2 < 0.001 
from the chosen profile. The left panels indicate to 
which models these profiles correspond. This is done 
by assigning to each model the same relative position 
as in the "model chart" shown in Fig. [31 The "model 
chart" is repeated four times: from left to right the 
models B_0, B_l, C_0, and C_l are represented (B: 
ballistic expansion, C: constant velocity, 0: no tur- 
bulence and 1: with turbulence ). The red square 
corresponds to the position of the chosen model, and 
the green squares correspond to all the models that 
give a profile that differs only by Di^2 < 0.001 from 
the chosen one. The result is spectacular. There 
are, in general, many models that have almost in- 
distinguishable profiles, and these models may have 
turbulent or non-turbulent velocity fields. If relying 
on just one line profile, the number of solutions may 
be very large. If, however, several lines are analyzed 
at the same time, the number of solutions is con- 
siderably reduced. If morphology and line profiles 
can be analyzed together, the number of solutions is 
further restricted. 

The main virtue of this exercise is to visualize 
that the possible solutions may actually be quite nu- 
merous. Because of the very coarse grid of models 
in the atlas, we obviously miss many other combi- 
nations of parameters. This shows that any solution 
obtained by model fitting cannot be regarded as nec- 
essarily correct. We emphasize that all of our mod- 
els have very simple density laws and geometries. If 
we were to consider more complex morphologies the 
number of solutions would grow accordingly. 

Notice that no instrumental profile has been ap- 
plied, so that the user can applyhis own instrumen- 
tal profile to themodeldata available from the server. 
Taking into account the instrumental profile would 
increase the number of profiles similar to a given one. 

It is straightforward, with C3D, to perform the 
same experiment for different velocity laws, and thus 
to visualize the different solutions corresponding to 
an observed line profile. 




Fig. 5. Left: In orange: position, in the "localization 
map", of the "original" model; in cyan: position, in the 
"localization map" of the models with line profile differ- 
ing from the original one by by D 1,2 < 0.001- Right: the 
corresponding line profiles. Each row corresponds to a 
different line, from top to bottom: H/3, He I 5876A, Hell 
4686A, [N II] 6584A and [O III] 5007A. 
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4. HOW GOOD ARE SIMPLE RECIPES TO 
ESTIMATE EXPANSION VELOCITIES? 

In many instances, either because there is no spa- 
tial information on the nebula or because it would be 
too time consuming to make an appropriate line pro- 
file fitting for a sample of nebulae taking into account 
all possible solutions, one relies on simple formulas 
to estimate the expansion velocity. 

Taking our atlas as a reasonable representation of 
the different possible cases in nature, we can estimate 
the errors and biases in the formulas commonly used 
to estimate expansion velocities. The distribution of 
morphologies of real PNs is obviously not compara- 
ble to the one we adopt here, but the global results 
of this work are still valid. 

To do this, we construct histograms of the 
HWHM velocities (expressed in units of mass- 
weighted mean velocities, vm), grouping the mod- 
els according to various criteria. Figure [6] upper left 
shows the histogram for all the virtual observations 
of the up line. In Fig.[6]upper middle, the virtual ob- 
servations are grouped into double peak profiles (red) 
and single peak ones (green). In Fig. [6] upper right, 
the models are grouped by morphology: spherical 
ones are in green, bipolar ones are in red, elliptical 
ones are in black. In Fig. [6] lower left, the group- 
ing is by apertures: red is for apertures covering the 
entire objects, green is for long slit apertures, blue 
is for small apertures, whether the aperture is cen- 
tered or off-center. Fig. [6] lower middle distinguishes 
objects with a various velocity fields. Finally, Fig. [6] 
lower right distinguishes the nebulae by orientation: 
face-on ones are in red, pole-on are in green, tilted 
are in blue. In the top of each panel, the crosses in- 
dicate the median value, and the extremities of the 
inclined segment indicate the quart iles, the color of 
each cross being the same as that of the correspond- 
ing histogram. 

It is striking that, in all the cases, the median 
values are displaced from 1 by less than 20% and 
the quartiles extend to 10 - 15% from the median. 
This indicates that such a simple estimate as HWHM 
is actually a quite accurate estimate of the mass- 
weighted expansion velocity. 

Figure[7lis analogous to Fig.[6l but now with half- 
width at 1/10 maximum (HWIOM) velocities instead 
of HWHM. While the histograms are again rather 
peaked, we note that the median value is significantly 
larger than 1 (typically 1.2 - 1.3). Therefore, the 
HWIOM velocity gives a more biased estimate of vm- 
Since Vhwiom is further plagued by the influence of 
noise to determine the position of the continuum, we 
conclude that Vhwhm is to be preferred to estimate 



the expansion velocity field. 

The parameters that have a large influence on the 
velocity determination are the ones that produce dif- 
ferent mean values. We find that the global shape of 
the profile (single or double peak) and the relative 
size of the aperture (small square, long slit, whole ob- 
ject) clearly affectes the inferred expansion velocity. 
On the other hand, for HWHM velocities, neither 
the real morphology of the nebula, nor the expan- 
sion velocity law, nor the orientation of the nebula 
have a strong influence on the velocity determina- 
tion. The latter parameters are those which are not 
known from the observer, while the first ones are 
always determined. These results attest the reliabil- 
ity of simple velocity estimator HWHM in the sense 
that invisible parameters do not bias the determined 
values. 

Table [2] lists the ratios of the velocities deter- 
mined by the HWHM and HWIOM methods rela- 
tive to the "true" values (mass- weighted velocities), 
in various cases depending on the shape of the pro- 
file (single or double peaked) or the relative size of 
the aperture used. The "errors" correspond to the 
25 and 75 percentiles. It is important to notice that 
this list of ratios is based on histograms constructed 
using the models from our atlas, where no attempt is 
made to reproduce the real distribution of PN mor- 
phologies and that those numbers are only indicative 
of the effects of various parameters but are not to be 
taken blindly as correction coefficients to apply to 
real observations. 

The effect of line splitting due to fine structure 
is obvious in the case of He I 5876A and He II 4686A 
lines, where the blueshifted low intensity component 
of the multiplet is still intense enough to be taken 
into account by the HWIOM method. The deter- 
mined expansion velocities are then determined 80% 
higher than the true values, while for unsplitted lines 
they are overestimated by 20 to 30%. 

The general conclusion is that expansion veloci- 
ties are reasonably well estimated from Vhwhm , es- 
pecially if one applies the correction given in Tab. [21 

5. SUMMARY 

Using the pseudo-3D photoionization code 
Cloudy _3D, we have constructed a grid of photoion- 
ization models of planetary nebulae with different 
morphologies, including elliptical and bipolar ge- 
ometries. We have then assumed different expan- 
sion laws, and computed line profiles correspond- 
ing to different orientation of the models consid- 
ered as virtual nebulae, and to different slit sizes 
and positions. The resulting atlas is available at 
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TABLE 2 

VELOCITY CORRECTION COEFFICIENTS. 



Line 


All 


Single Peak 


Double Peak 


Full object 


Long Slit 


Small Aperture 


HWHM 


H/3 


^'^^ -0.13 


0.84 


1 19 +0.15 
-^•-^^ -0.10 


r| 77 +0.06 

^' ' ' -0.03 


'-'•^^ -0.09 


1 20 +°-°^ 
^•"^^ -0.12 


He I 5876A 


J-.UO _o.i4 


-1 no +0-09 
J-.uz _o.i6 


1 +0.18 
-"-•-■-^ -0.13 


0.86 tUl 


1 0^ +°-°'^ 


1 SO ^^-^^ 

J-.OU _o.i4 


He II 4686A 


-0.09 


1.02 iro? 


-, 9-, +0.13 
^•"^^ -0.09 


0.96 iro^ 


1 10 +°-^'^ 

-0.07 


1 SO +°-^^ 

-0.14 


[Nil] 6584A 


1.05 t^;?? 


0.88 in^ 


-'-•'-'^ -0.05 


0.88 iro^ 


1.05 t^;^o? 


-1 -19 +0.06 
^'^^ -0.06 


[0 III] 5007A 


0.98 t'o-'ol 


0.85 ^Q-Jg 


-'-•'-'^ -0.09 


0.86 iro? 


0.98 troi 


1 16 +°-°^ 
-'-•-'-'^ -0.10 


HWIOM 


Rp 


i.^o _o.08 


-, 07 +0.08 
-0.09 


-0.08 


-1 90 +0.09 
^'^^ -0.02 


1.43 l^;^o's 


1 60 

i.OU _o.09 


He I 5876A 


i.OO _o.i2 


-1 ro +0.09 

-'-•^'^ -0.10 


-1 ro +0.23 
1.05 _o.i3 


i.^O _o.09 


1.58 1^;?^ 


-1 7r +0.07 
-■-.'^ -0.15 


He II 4686A 


-1 o-i +0.18 
J-'^J- -0.11 


1 7Q +0-19 
-0.17 


-I nA +0.14 


^'^^ -0.03 


1 7Q +0-13 
' ^ -0.08 


-I QO +0.14 


[Nil] 6584A 


-, 90 +0.14 
-'-•^'^ -0.09 


-, 90 +0.07 
i.^O -0.13 


-1 -.0 +0.20 


-, 17+0.11 

-0.12 


-1 97 +0.10 
-0.13 


-1 01 +0.26 
^'•^^ -0.17 


[0 III] 5007A 


-1 90 +0.10 
-0.13 


-1 90 +0.09 
_o.i2 


-1 90 +0.12 
_o.i3 


-1 9-1 +0.07 
-0.16 


1 SO +°-^^ 

-0.15 


1 40 +°-^^ 

-0.26 




R: Double peack 



Shape: R: Bip, G: Sph, B: Ell 





0.0 0.5 1.0 1.5 2.0 



0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



Aperture: R: Full, G: Slit, B: Squared 



Exp Vel: R:B_0,G:B_1 ,B:C_0,Blac:G_1 Orientation: R: Face-on, G: Pole-on, B: 45deg 





0.0 0.5 1.0 1.5 2.0 



0.0 0.5 1.0 1.5 2.0 



Fig. 6. Histograms for the HWHM-velocity determination normalized to the mass-weighted real values, for 
See text for a detailed description 



the line. 
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R: Double peack 



Shape: R: Bip, G: Sph, B: Ell 
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0.0 0.5 1.0 1.5 2.0 



0.0 0.5 1.0 1.5 2.0 



Aperture: R: Full, G: Slit, B: Squared 




Exp Vel: R:B_0,G:B_1 ,B:C_0,Blac:G_1 Orientation: R: Face-on, G: Pole-on, B: 45deg 
1.0 r 
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/ 
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0.0 0.5 1.0 1.5 2.0 



Fig. 7. As Fig. [6] but for HWIOM- velocities. 



http://132.248. 1. 102/Atlas-profiles/ . In its 

present version, it is already useful to visualize the 
variety of profiles that can be obtained according to 
the physical properties of the nebulae and to the ob- 
serving conditions, and to help in the interpretation 
of observational data. 

We have shown that in many circumstances, line 
profiles are degenerate, and the recovery of the true 
geometry and velocity field from observations re- 
quires lines from ions with different masses and dif- 
ferent ionization potentials. 

We have also shown that the empirical way to 
measure mass- weighted expansion velocities from the 
HWHM of observed line profiles is accurate within 
about 20% in the range of geometries and velocity 
fields considered in the present version of the atlas. 
In the case of observations of distant nebulae which 
are entirely covered by the slit, the unknown geom- 
etry and orientation of the nebula do not alter the 
measured velocities, when considered statistically. 

We plan to integrate our atlas in a virtual ob- 
servatory environment. For the time being, the 
Cloudy _3D suite and the associated VISNEB tool 
are available on request from CM. and may be used 
as a help to interpret observed expansion profiles. 

One application of the atlas will be to compare 
the synthetic line profiles and PV-diagrams with ob- 
served data such as the SPM Kinematic Catalogue 
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